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We propose a non-perturbative attempt to solve the kinematic equations for spiral waves in 
excitable media. From the eikonal equation for the wave front we derive an implicit analytical 
relation between rotation frequency Q. and core radius Rq. This Q.{Rq) dependence improves the 
result by Keener and Tyson for spiral waves pinned to a circular defect with Neumann boundaries 
at the periphery. Simultaneously, for free, rigidly rotating spiral waves our analytical prediction is 
in good agreement with numerical solutions of the linear eikonal equation not only for very large 
but also for intermediate and small values of the core radius. We discuss the reasons why the ansatz 
fails to correctly describe the result for the dependence of the rotation frequency on the excitability 
of the medium as obtained from the numerical solution of the eikonal equation for the wave back. 
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I. INTRODUCTION 

Spiral waves are a common occurrence in excitable 
media. They have been observed in systems as diverse as 
catalytic surface reactions |T], the Belousov-Zhabotinsky 
chemical reactions |21 13] and social amoeba colonies 
They play an important role as pathological states of 
action potential propagation in cardiac tissue and are 
thought to be the precursor of ventricular fibrillation [5J . 

In the simplest free spiral rotates rigidly with 

a frequency uj while its tip describes a circular trajectory 
with radius tq, called the spiral core radius. From exper- 
iments and numerical simulations it is well-known that 
spiral waves select their own, unique asymptotic wave 
shape and rotation frequency. Thus, independently on 
the method of initiation, coexisting (and non-interacting) 
spiral waves in a spatially uniform excitable or oscilla- 
tory medium exhibit the same wave length, core radius 
and rotation frequency after all transients have died out. 
Exceptions to this rule are known only for media with 
anomalous dispersion of periodic wave trains . 

A theoretical description aims to understand the un- 
derlying selection principle determining shape, rotation 
frequency and core radius of spiral waves. One successful 
approach is the free boundary or kinematic approach 
[71 |S] which reduces the nonlinear reaction-diffusion 
equations to simpler equations describing the motion 
of interfaces separating excited from resting states. In 
contrast to earlier works, which mapped wave front and 
wave back onto each other [71151 [10]; it became clear that 
one has to solve equations for both the front and back 
interface of a spiral to yield a self-consistent solution for 
LJ and tq. 

Within the free-boundary approach, the pattern selec- 
tion problem for wave segments which are unsta- 
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ble solutions of the reaction-diffusion system, has been 
solved. These patterns undergo translational motion in 
an unbounded medium. The properties of the medium 
were expressed by a single dimensionless parameter B 
that can be interpreted as a measure of the local ex- 
citation threshold which increases with B while the ex- 
citability decreases. Wave segments exist in a finite range 
< -B < Be of B-values. 

In the limit i? — >■ 0, wave segments correspond to 
motionless circular spots [12J. At the upper boundary 
of the existence range, for B — > i?c~0.535, they merge 
with spiral waves in a separatrix solution known as the 
critical finger |13] . The critical finger is an half-infinite 
plane pulse with an open end that can be regarded as 
a spiral wave rotating with zero frequency around an 
infinitely large core. 

Spiral waves with B < Be close to the critical finger 
and their transition into meandering spiral waves were 
analytically investigated by Hakim and Karma [2] 
applying methods of singular perturbation theory. For 
B > Be, only retracting fingers were found to exist 
because the excitability of the system is too weak to 
support wave segments and spiral waves. The entire 
range -Bmin l£ B < B^ for which spiral waves exist 
as a solution to the kinematic equations was studied 
numerically by Zykov in [HUTS]. For B B,nin ~ 0.211 
the spiral core radius tq decreases to zero and for 
B < Bmin, rigidly rotating spiral waves cease to exist. 

Solving the kinematic equations in a disk with a Neu- 
mann boundary led to the discovery of boundary spots. 
Boundary spots are unstable wave solutions to reaction- 
diffusion equations that rotate at a lower frequency than 
spiral waves. Furthermore, in contrast to spiral waves, 
boundary spots do not extend over the entire disk but 
remain spatially localized close to the disc boundary [17J. 

In this work, we follow an analytical non-perturbative 
approach that goes back to a classical paper by 
Burton, Cabrera and Frank [TS]. This approach is 
non-perturbative in the sense that it does not rely on 
a small parameter for a perturbation expansion. These 
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authors used an ansatz function for the wave shape 
to calculate the rotation frequency of spirals rotating 
around a core with zero radius. Keener [12] and Tyson 
and Keener fH^ extended this ansatz to spirals pinned 
to a circular hole with a no-flux boundary condition 
at the periphery. In this case the selection problem 
turns out to be much simpler than for free spiral waves 
because the rotation frequency can be determined from 
the equation for the front interface alone while the core 
radius is given by the radius of the Neumann hole. 

Below, we present a non-perturbative approach which 
does not only improve the result obtained by Keener 
and Tyson for pinned spiral waves, but also works quite 
well for free spirals. In Sec. |TT] we state the equations 
of the free-boundary approach and review existing solu- 



tions. Our ansatz is introduced in Sec. Ill The analytical 
results for free and pinned spirals are presented in Sec. 
|IV| and Sec. [V] respectively, and compared to numeri- 
cal solutions of the kinematic equations. We end with 
discussion of the results, conclusion and outlook in Sec. 
lYll 



II. KINEMATIC EQUATIONS 

We consider a standard activator (u) -inhibitor (v) 
reaction-diffusion systems of the form 



dtu = e\'^u + f{u,v) /e, 
dtV = g{u,v) , 



(1) 
(2) 



where the dimensionless parameter < e ^ 1 is a mea- 
sure for the time scale separation between activator and 
inhibitor and serves as a small parameter for a pertur- 
bation expansion. We neglect inhibitor diffusion and 
scale space accordingly so that the activator diffusion 
coefficient is equal to e. The u nuUcline obtained from 
/ (m, d) = is assumed to be S'-shaped in the (u, v) plane. 
A simple choice for the functions / and g is given by the 
FitzHugh-Nagumo kinetics 



g {u,v) = u — 



(3) 
(4) 



with a unique, linearly stable rest state uq = S, vq — 
35-6^. 

If e is small, a traveling pulse can be regarded as con- 
sisting of two separate spatial regions: an excited region 
where the value of the activator is large and the 
inhibitor is rising, and a refractory region where 
the activator value is small and the inhibitor is decaying. 
This behavior is described by the outer equations Eqs. 
([ij , Q , which in lowest order to e read j8] 

= f{u^iv),v), (5) 
dtV = g (li* (u) , u) inV^. (6) 

Here, u+ (v) and (v) denote the largest respectively 
smallest root of f {u,v) — which the activator follows 



in the excited respectively refractory region. The two 
regions I?+ and 'D~ are separated by a front and 
a back (— ) interface, where the activator value changes 
very fast from a low to a high value and the other way 
round, respectively. These interfaces can be regarded as 
fronts traveling with velocities c*. They are solutions 
to the inner equations, obtained from Eqs. ([ij, ^ 
by a change of scale in time and space proportional 
to e. The expression for the front velocity together 
with Eq. ([6 1 and appropriate periodic boundary 
conditions yield the dispersion relation for a periodic 
pulse train, i.e., the dependence of the propagation ve- 
locity c on the period length L to lowest order in e [8, 19J. 

In two spatial dimensions, the shape of the front 
(+) and back (— ) interfaces for rigidly rotating spiral 
waves are conveniently parametrized by 6^ (r) using 
polar coordinates 



x± (r,<) 
2/± {r,t) 



r cos {9^ (r) — uit) 
r sin {9^ (r) — ujt) 



(7) 



In Eq. Q, w > is the rotation frequency of a spiral 
wave rotating counterclockwise. The inner equations in 
two spatial dimensions provide a relation between the 
normal velocity of the moving front and back interface 
and its local curvature [12] , the so-called linear eikonal 
equation 



(8) 



Here denote the inhibitor level at the interface, and 
(v^) is the velocity of a planar front moving through 



a medium with a constant inhibitor value 



Similar 



as for a one-dimensional pulse train, Eq. (|6| yields to- 
gether with the condition of periodicity in 9 an expression 
for {v^) ■ This constitutes the so-called wave front in- 
teraction model. The interaction between wave front and 
wave back is mediated through the dependence of i;+ and 
w~ on the positions of both front and back interface. 
With the chosen parametrization, the curvature is 
expressed as 



K± (r) 



9^' (r) 



(d/dr) (r6l±' (r)) 



(l + M±'(r))')'^' (l + (r0±'(r))') 



3/2 ' 



(9) 



and the normal velocity is given by 



(l + (r0±'(r))') 



1/2- 



(10) 



Eq. Q has to be supplemented with appropriate bound- 
ary conditions. For a rigidly rotating free spiral wave, 
front and back interface meet continuously at the apex 
r — ro of the spiral, i.e.. 



(r-o) = 0, 



(11) 



3 



where we fixed an arbitrary initial phase of the spiral to 
be zero. The apex is the point of closest approach of both 
interfaces to the center of rotation (compare Fig. [l]). At 
the apex, the normal velocity is zero, (rp) = 0. Both 
interfaces approach the apex tangentially to a circle with 
core radius tq, so that 



5+' 



(ro) = -e-' (ro) = (X3. 



(12) 



This circle is considered as the spiral core with ro being 
the core radius. 

Far from the core, front and back interface behave as an 
Archimedean spiral. 



(r) ^ r, (r — )■ oo) . 



(13) 



Eqs. (Ill, (12), (13) fix six boundary conditions for 



two coupled second order ordinary differential equations 
(ODEs) Eq. ([S]). Four boundary conditions are neces- 
sary to determine the four integration constants of these 
ODEs. The remaining two are used to determine two 
unknown nonlinear eigenvalues introduced as parameters 
in the eikonal equation and the boundary conditions: the 
rotation frequency lu and the spiral core radius tq. The 
full wave front interaction model, as given by Eq. ^ 
together with the linear eikonal equation Eq. ([s]) was 
solved numerically by Pelce and Sun in |20 without any 
further approximations for a piecewise linear activator 
kinetics. 

Assuming that the inhibitor value v stays always close to 
the stall level v = Vs given by (vs) — 0, Eq. ^ can 
be simplified [2] 



dtv 
dtv 



Wo 



in V 



with the abbreviations 
1 



TR = 



duf 



dugdvf ~ d^gduf 



(14) 

(15) 
(16) 

(17) 
(18) 



U — U (Vs), V—Vs 



This approximation assumes a linear rise of the inhibitor 
during the excited period on a time scale of the order , 
followed by an exponential decay during the refractory 
period on the time scale tr. 

Spiral waves close to the critical finger have a diverg- 
ing refractory period, so that the inhibitor value has 
already decayed to its rest value, = vq, everywhere 
along the front interface. In this case depends linearly 
on the angular pulse width A6 (r) = 6^ (r) — 6~ (r), and 
the expressions for become particularly simple 



{v^) = c, 



c+-{e+ (r) 

UJ 



^)) 



(19) 
(20) 



c = a {vs — Vq) > corresponds to the velocity of a front 
solution of the inner equations moving through a medium 
with a constant inhibitor at its rest state v — Vq. Note 
that the eikonal equation for the front interface decouples 
from the equation for the back, while the back interface 
interacts with the front interface via a term that is linear 
in the pulse width. The single kinetic parameter b = 
a/Te > is a measure for the strength of this interaction. 
For FitzHugh-Nagumo kinetics according to Eqs. (|3|, 



((4|, we find a = 1/a/2, Vg = 0, Te 



and Tr = 6. 



V3~S 

Hakim and Karma |14j used singular perturbation theory 
to simplify the eikonal equation Eq. ([8|. Taking into 
account Eqs. ([T9|, (pOl follows 



Cn ir)^c- en+ (r) , 

c" (r) = - (61+ (r) - e- (r)) - c - e/c" (r) . 



(21) 
(22) 



Again, e serves as the small parameter for a singular per- 
turbation expansion. For both the front and back inter- 
face three scaling regions were identified: the spiral tip 
region near to the core, an intermediate region, and one 
region sufficiently far from the core where curvature ef- 
fects are less important. The outer equations valid in the 
far core region are Eq. ([s]) with e = 0. Its solution 



(r) = r6'+ ' (r) 

T^mv \ / inv V / 

er (r) = e+ (r) - 

inv V / ^inv V / 



2CL0 



(23) 
(24) 



describes the involute of a circle of radius ro which 
asymptotically transforms into an Archimedean spiral for 
r oo. Imposing the boundary condition 



9+ ' 

mv 



(ro) = 0, 



gives 



(25) 



(26) 



while the spiral core radius ro is left undetermined. The 
boundary condition 



9± 

mv 



(ro) = 



(27) 



can only be fulfilled either at the front or the back inter- 
face but at not at both interfaces together. 
The behavior in the tip region is described by the equa- 
tions for the critical finger. While the equation for the 
front can be solved analytically [111, no analytical so- 
lution is known for the back. Matching the analytical 
solutions for the front interface in the tip, intermediate 
and far core regions. Hakim and Karma succeeded to 
derive analytically a relationship between rotation fre- 
quency UJ and core radius ro. Matching the solutions for 
the back interface and using stability arguments, an ex- 
pression for the core radius ro involving two numerically 
determined constants was obtained. Together, these two 
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relations yield the desired dependence of the spiral wave 
frequency w on the kinetic parameter b. It should be 
emphasized that these solutions are only valid for small 
e and for spiral waves close to the critical finger which 
have a diverging core radius. 

The eikonal equations for wave front and back, Eqs 
( 22 I . can be rescaled by introducing dimensionless quan- 



(211 



titles according to 



Re 

c 



and rescaled shape functions 9* as 

'Re 



(r) = 



= 0" (R) 



cCt, (28) 



(29) 



and 



*± (R) = i?e±' (i?) 



(30) 



In the following all rescaled dimensionless quantities will 
be denoted by upper case letters, while lower case let- 
ters are used for dimensional quantities. In dimensionless 
terms Eq. (|2ll) and Eq. (p2| read 



C+ (i?) = 1 - K+ {R) 
B 



(31) 



CniR) 



-1 



e+ (R) - e- (i?)) - (R) . (32) 



Note that the small parameter e as well as the propaga- 
tion velocity c in the eikonal equations have been elimi- 
nated under the rescaling. Applying a numerical shoot- 
ing method, Zykov [151 ITB] demonstrated the existence of 
spiral wave solutions to these equations in a certain inter- 
val -Bmin ~ 0.211 < B < Be ~ 0.535 of the dimensionless 
excitability parameter B and determined the desired uni- 
versal relationship Q. (B). 

At B = i?min ~ 0.211 the shape of the front interface is 
identical to that obtained by Burton, Cabrera and Frank 
(BCF) [18, for a spiral wave with zero core radius rotat- 
ing at frequency ft « 0.331. The back interface results 
from turning the front interface clockwise around an an- 
gle AG (i?) = 6+ {R) ~ Q- (R) = TT. In the other hmit, 
for B<Bc~ 0.535, the analytical results of Hakim and 
Karma for spirals with diverging core radius are recov- 
ered [16]. 

The front interface of spiral waves with finite core ra- 
dius i?o interacts with the back interface of the wave 
ahead because it does not propagate into a fully recovered 
medium. In general, the front inhibitor level u"*" depends 
on the radial coordinate r. However, assuming a constant 
value of the inhibitor given by the dispersion relation of a 
one-dimensional periodic pulse train and using a slightly 
different rescaling, the dimensionless eikonal equations 
Eqs. (|3l]), (|32| can still be used [T5] . 

This approximation clearly separates the two physical 
mechanisms underlying the frequency selection for spiral 
waves: 

I. The interaction of the front interface with the back 
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Figure 1; Close-up of the spiral tip region. Front (£»+ (r), 
solid line) and back [9~ (r), dashed line) interface separate 
the shaded excited region from the white refractory region 
O". At the core radius ro, the point r+ at the wave front and 
the inflection point ri at the wave back, approximate and 
exact analytical solution of the free-boundary problem have 
the same leading order asymptotic expansions. 



interface of the preceding wave essentially leads a front 
moving through a only partially recovered medium. This 
in turn leads to a slower velocity of the front as ap- 
proximately given by the dispersion relation of a one- 
dimensional periodic pulse train. 

II. The interaction between of the back interface with the 
front interface within the same wave is proportional to 
the angular pulse width AO (R) = 9+ (i?) - 9" (R) and 
characterized in strength by the dimensionless parameter 
B. 

That the kinetic characteristics of the medium can be 
lumped together into a single parameter B simplifies the 
analysis of the existence of spiral waves significantly. 

Front and back interface of spiral waves pinned to a 
circular Neumann hole of radius r+ approach the hole in 
radial direction, so that 6^ fulfills the boundary condition 



Q±' 



(33) 



which implies that the spiral arm is orthogonal to the 
hole. Note that a corresponding point r — can always 
be found at the front interface of a freely rotating spiral 
wave, see Fig. [T] The kinematic equation for the front 
interface of a pinned spiral wave was studied analytically 
by Keener and Tyson |H1 HH] • These authors determined 
the asymptotic behavior of solutions to Eq. (21 1 for the 
front interface as r — >■ r_|_ and r — >■ cx). An ansatz showing 
the same asymptotic behavior and involving several con- 
stants was used. Comparing the asymptotics of ansatz 
and solution, they were able to determine the constants 
of the ansatz and finally derived a relation between rota- 
tion frequency a; and hole radius r_|_. 
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III. ASYMPTOTIC SOLUTIONS AND ANSATZ 0.5 



In this chapter we obtain asymptotic sohitions to the 
hnear eikonal equations Eqs. (21 1, (22 1 for the front and 



back interface of rigidly rotating spiral waves valid at the 
spiral core, r — ro, at the point r = r+and very far from 
the core as r — >■ oo. Additionally, the existence of an 
inflection point at r = ri at the back interface is taken 
into account. Then we present an ansatz for the interface 
shape 6^ (r) that reproduces in leading order all these 
asymptotic expansions correctly. See e. g. [21 how to 
compute asymptotic solutions to differential equations. 



A. Asymptotic solutions to the linear eikonal 
equation 

Far from the core the shape of the interfaces is asymp- 
totically Archimedean, i.e., 



(r) = r + O (log (r)) , r — J' oo. 



(34) 



At the spiral core, the asymptotic behavior that fulfills 
the two boundary conditions Eqs. (Ill, ( 12 1 is given by 
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Figure 2: Dimensionless rotation frequency S7 over dimension- 
less spiral core radius Rq for rigidly rotating spiral waves. The 
analytically obtained result according to Eq. ( 49 1 (black solid 
line) is compare d to num erical solutions of the linear eikonal 
equations Eqs. (311, (321 (red dashed line). The blue dotted 



line shows the analytical result for diverging core radius [14] . 

B. Asymptotics of the ansatz 

We use the ansatz 



(r) 



± 



--\Jt - J'o + O (f - ''o) , T^Tii. 



(35) 



Finally, at the distance r = r+ (compare Fig. [T]), an 
asymptotic expansion for the front interface is available 
which reads [5J [H] 



r+e 



(r - r+) O [(r - r+)^^ , r ^ r+. 

(36) 



The back interface of spiral waves always exhibits an in- 
flection point at r = ri where the curvature vanishes 



V-^ns (j) = vQtJ (r) = A± . (40) 

which yields for the interface shape 9^^^ (r) 



.± (,) = / df 



A± ( \/ — r?: — arccos 



(7))- (") 



Our ansatz produces the correct leading order asymp- 
totics for r— > 00 



K- (n) = 0. 



(37) 



In polar coordinates the inflection point is not easy vis- 
ible, see Fig. ^ while it appears clearly pronounced in 
Cartesian coordinates used in Fig. |4] From the eikonal 
equation for the back Eq. ( 22 I follows at the inflection 
point 



c„ (ri) = c (ri) 



(38) 



which using Eq. ([10| leads to the following expression 
for the parameter b 



b = 



riUJ 



(39) 



9± 

ans 



(r) = A±r + O (1) , r ^ 00, 



(42) 



and for r tq 



+ O ((r - ro)'/') , r ^ ro. (43) 

At the point r = r_|-, the ansatz for the front interface 
displays the asymptotic behavior corresponding to a Neu- 
mann boundary 



9+ 

ans 



(r) = (r -r+)+0 ((r - r+ f) , 



(44) 
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Figure 3: Rotation frequency versus excitability param- 
eter B in dimensionless units. The red dashed and dash- 
dotted lines represent the branches of spiral waves respec- 
tively wave segments obtained numerically from the linear 
eikonal equations Eqs. (311,(321. The two branches merge for 
B = Be ~ 0.535 (critical finger |13p. The analytical approx- 
imation is shown by the black solid line where the thick seg- 
ment corresponds to the frequencies below Q — Qbcf ~ 0.331 
(Burton-Cabrera- Frank limit) . The blue dotted line shows the 
analytical result from |14) . 



For r_ < Tg, 9^^^ exhibits an inflection point at a point 
ri > ro- Thus, valid solutions for the front and back of 
a free spiral can only be found if 



> ro > r_ > 0. 



(45) 



In the following we consider only solutions which fulfill 



condition Eq. ( 45 1 



Equating the six asymptotic expansions of the ansatz 
with the six asymptotic solutions to the linear eikonal 
equation, we are able to determine four unknown param- 
eters A± and r±, as well as the relations uj over ro and 
To over b. 



IV. RESULTS FOR FREE SPIRALS 



A. f7 over Ro 



Comparing Eq. ^ with Eq. (|42|), and Eq. (|34| with 
Eq. (|43|, we get 



r±^ \ ro{ro± 



y cro + e 



(46) 
(47) 



Equating the asymptotic expressions for r ^ given by 



Eqs. ( 44 1 and ( 36 ) we obtain an implicit relation between 



frequency uj and core radius ro 



cro 



cro+e 



J -'V^'(cro + e) ^ (48) 



c^ro 



In rescaled form, Eq. (48 I reads 



in 



i + RoVWTin 

from which we find 
1 1 



-2Jn3^/Ro + l. (49) 



2V2 



'R 



^ 



and 



n = 



1 



^0 2K 



3/2 



0{R, 



-7/4 



, ^0 



(50) 



(51) 



Thus, we obtain a finite rotation frequency for vanishing 
core radius. However, with 17 (i?o — 0) — 1/2 we miss 
the Burton-Cabrera- Frank limit Vl « 0.331. 



Eq. ( 49 I can be written as 



+ {Rfi + 2Vn}j - i?2 (1 + i?o) = (52) 
where we have we have introduced the abbreviation 

(53) 



which can be used to rewrite R± from Eq. (47) in the 
form 



R± 



Ra\ Rf)^ 



1 



VI + Ra^ 



Ro\ll±^. (54) 



From Descartes rule of signs we conclude that the 
number of positive real roots of the sextic polynomial 



Eq. (52) is one. Thus, ft over i?o has only one physically 



meaningful branch. 

Our theoretical prediction for D, (i?o) matches well with 
direct numerical solutions of the rescaled eikonal equa- 
tions Eqs. (31), (32) (compare black solid respectively 
red dashed line Fig. |2| where in particular the good 
agreement for intermediate and even quite small core 
radii is remarkable. 



Perturbatively, we can derive from Eq. (48 ) an explicit 
expression for w(ro) valid for small e, 



w (ro) 



„l/2 



2r, 



3/2 



^1/2+0(^,3/4^ 



(55) 



Hakim and Karma found by singular perturbation theory 
for very large core radius Rq [TJ 



^HK (ro) = - + 21/3^,2/3 



'"0 



(56) 
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Upon rescaling (39 1 and comparing it with the ansatz, 



Figure 4: Front and back interface of free spiral waves in 
the core region. Black solid and dashed lines show analytical 
approximations for wave front respectively back as obtained 
from Eq. ||69| using Q.{Rq) given by Eq. (|49l. Red dotted 



and dot-dashed lines are plots of the corresponding numerical 
results. The rotation frequency was fixed to a) f2 = 0.137, b) 
Vl = 0.102, c) n = 0.219. The analytically obtained core radii 
(black circles) are a) i?o = 3, b) iio = 5, c) _Ro = 1- 



which in rescaled form reads 

1 



^^HK (-Rq) 



r)l/3 '^l 

^ 5/3 ' 





(57) 



where a\ = —2.3381 denotes the first zero of the Airy 
function Ai{x). The dependence Eq. (57 1 corresponds 
to the blue dotted hne in Fig. [2]Note that for e = 0, we 
find 



ro 

r+ = ro, 

i.e., our ansatz reduces to the involute spiral 



c->0 




(58) 
(59) 



(60) 



Eq. (60) gives the correct solution of the linear eikonal 



equation for e = 0, compare Eq. (23), which motivates 
our ansatz a posteriori. 



n over B 



we can express B in the form 



B 



^fl,f(e-„3'(i?i))^+i 



©ins (^l) — 6ans (^l) 



Here, Ri is given by 



K- (Ri) = 



or explicitly 



Rfe-J iRiY + RiQ-J' (i?i) + 2e 



(61) 
(62) 

(i?i)-0. (63) 



If we numerically solve the algebraic Eqs. (|6T|) and (63) 



together with the ansatz Eq. (|41j) and use the relation 
for 51 over Rq Eq. ( [49| , we obtain the dependence fl (B) 
plotted as black solid line in Fig. [3j Analytically we 
can determine _B as a function of O in the limit il — >■ | 
corresponding to i?o — >■ 

lim B (Q) 



1 



0.437171, 



a^i ' • 4arcsec (21/4) 
and in the limit fl — )• corresponding to Rq 
lim B (fl) =0. 



(64) 



(65) 



From the obtained il [B] dependence we conclude that 
spiral wave solutions to the linear eikonal equation exist 
only within a finite range of B values. The red dashed 
line in Fig. [3]shows f2 {B) calculated numerically for free 
spiral waves in [T5|. The red dot-dashed line in Fig. [s] 
we have added for completeness. It describes the branch 
of unstable wave segments. At i? = i?c ~ 0.535 The 
two red branches merge in the critical finger. The blue 
dotted line in Fig. [3] shows the analytical result derived 
by Hakim and Karma |T3] 



flHK [B) = 



V2 



B -B, 
Kai 



3/2 



(66) 



In 



which is valid close to the critical finger B < B, 
Eq. (41 1 K K 0.630 and Be are numerically determined 



constants. 

For the interpretation of the obtained ft (B) dependence 
one has to bear in mind that the linear eikonal equation 
has rotating wave solutions of different type. In a disk of 
finite radius Rd with Neumann boundary conditions at 
the periphery, stable spiral wave solutions coexist with 
unstable boundary spots for a given value of B. The 
latter represent localized unstable solutions attached to 
the boundary that rotate at lower frequency as compared 
to the corresponding spiral wave |15l [Tf]. Increasing B, 
spiral wave branch and boundary spot branch coalesce 
in a saddle-node bifurcation at a certain value Bar (Rd) 
beyond which both spiral waves and boundary spots 
cease to exist. If the disc radius goes to infinity, the 
upper, high-frequency Q (B) branch describes free spiral 



8 



waves while the lower ft (B) branch corresponds to wave spiral wave 

segments which can be imagined as boundary spots 

"attached to infinity". Upper and lower branch meet q± ^jj") 

in a saddle-node bifurcation forming the critical finger 

which gives the upper bound for the existence of spiral 

wave solutions in the unbounded medium. 



i?2 - R^n 



Rnil ± . arccos — 



(69) 



In Fig. ([3]) the analytical approximation for fl (B) (black 
solid line) intersects the numerically calculated spiral 
wave branch (red dotted line). Imaginable would be 
that within an improved ansatz this intersection point 
shifts down along the numerically calculated spiral wave 
branch towards the critical finger, and simultaneously 
the branches of the black curve converge gradually to 
their numerically computed counterparts. In summary, 
although the deviation between the analytically and the 
numerically obtained fl (B) is definitely large, the former 
appears to be comprehensible in the suggested interpre- 
tation. 

Analyzing the reasons for the discrepancy between the 
analytically and numerically determined n (B) depen- 
dence it is important to note that our ansatz fails to de- 
scribe solutions with small rotation frequency. The shape 
of the wave front 0+ (R) obtained from the rescaled lin- 
ear eikonal equation Eq. ( |31[ ) behaves in the limit i? — > 00 
as 



e+ (i?) = -nR-n\n{R) 



1 

m 



1 



i? -> 00. 



(67) 



As B ^ Be and SI — > 0, all terms linear in 17 vanish while 
the term ~ l/f2 grows indefinitely. Therefore, expansion 
Eq. (67 1 breaks down close to the critical finger. The 



correct leading order asymptotics for the critical finger 
reads 



9+ (i?) = 2 



In(fl) 
R 



ln(i?) 
R 



, R ^ 00, B = Be 



(68) 



where 17 as a function of core radius Rq is given by Eq. 

4] we compare the analytical prediction 
utions of the dimensionless eikonal equa- 
"32I for three given values of the ro- 



(491. In Fig. 



to numerical sol' 
tions Eqs. ( [3l| , 
tation frequency 17 in order to avoid the inaccuracy in 
the analytical relation 17(5). Fig. |4] shows good over- 
all agreement between theoretical and numerical results. 
In particular, the front interface is nicely described by 
the ansatz although the analytical approximation always 
slightly overestimates the core size. However, the analyt- 
ically predicted pulse width 



arccos 



Rq 



(70) 



displays deviations that increase for large core radii. In 
fact, asymptotically we find for large R 

Ae.„.(i?) = ^+o(^),i?-oo, (71) 

while the asymptotic solution to the eikonal equation 
yields 

Ae (R) = e+ (i?) - e- (r) 

^l + Ks)-'*""- '''' 

If we plug the asymptotic expansion for the rotation fre- 
quency 



17 



1 

i?0 



1 
i?0 



(73) 



The reason for the breakdown of the expansion Eq. 
(67 1 for 17 — > is that polar coordinates are a conve- 



into (72 



nient parametrization for spiral waves but a bad choice 
for wave segments and critical finger which are better 
parametrized in Cartesian coordinates. Note, that the 
asymptotics for r — >■ tq, Eq. (351, remains valid for the 
critical finger and wave segments. In other words, the 
tip region of the critical finger and of wave segments is 
correctly represented by our ansatz, however, it fails in 
correctly predicting the whole shape of the front interface 
for wave segments and the critical finger. 



) we get a different asymptotic behavior for large 
ii as compared to (71 1. This difference explains 
the decreasing agreement between analytically and 
numerically calculated back interface with increasing 
Rq. Moreover, because B is intimately connected with 
the pulse width, it is another reason for the failure of 
the analytical 17 over B relation in the vicinity of the 
critical finger. 

The opposite limit, Rq — 0, implies 17 — 1/2, B — > 
1/ (4arcsec (2^/^)) « 0.437171 for the ansatz solution. 
This yields an Archimedean spiral 



C. Analytical approximation for the spiral shape 



lim (i?) 



-2(^T^) 



(74) 



Our ansatz leads to the following analytical approxima- 
tion for the front and back interface of a rigidly rotating 



with a back interface identical in shape to the front inter- 
face but turned by an angle AGans {R) — tt. Numerically 



9 



solving the kinematic equations for zero core radius, i.e., 
for w 0.331 and B = iJmin ~ 0.211, leads to a very 
similar result, compare red lines in Fig. [Sjand [TBI [^ . 
To find even better analytical estimates for B and ft in 
the limit Rq — > 0, we choose an Archimedean spiral with 
the correct leading order asymptotics as i? — >■ oo, i.e. 



(75) 



for the front interface, and assume an identical back in- 

terface turned by an angle A0ans {R) = tt = — [22\ . 

B 

With this ansatz, we minimize the functional 



S+= I dr {1 - K+ (r) ^ C+ ir)y 

4 - TT 





16""+^ 



log(4)-l 



(76) 



with respect to il. We obtain 



19 



(4 - tt) TT - 1 - log (4) w 0. 144 (77) 



for 



8 (4 - tt) 
197r 



0.339. 



The corresponding value for B is 



B 



2n 

TT 



2(4- 



197r3 



0.216. 



(78) 



(79) 



These values, though only approximately valid, display 
a relative error of less than 3% when compared with the 
numerical results of the BCF limit. We compare this 
Archimedean approximation for the spiral shape with nu- 
merical solutions of the rescaled eikonal equations Eqs. 
([31}, ([32]) in Fig. [5] 



RESULTS FOR PINNED SPIRALS 



Now, we consider a wave rotating around a hole of radius 
r = r+ with no flux boundaries at the hole boundary. 
Then, we are looking for solutions of Eqs. 
subjected to the following conditions 



(21), (22) 



e+ (r+) = 9+ 
0±'(r+)=O, 
0± (r) 



r, r oo. 



(80) 
(81) 
(82) 



Four out of these five conditions are necessary to deter- 
mine four integration constants of the ODEs. The fifth 
condition yields a relation between lu and r-|_. 



15 



10 



Y 



-10 
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Figure 5: Spiral shape in the BCF limit of zero spiral core 
radius. The back interface is identical to the front interface 
turned by an angle AQ = tt. Numerical solution of the eikonal 
equation (red) and Archimedean approximation given by Eq. 
( 75 1 (black) are plotted in such a way that they coincide far 
from the core. This leads to small deviations near to the 
center of rotation. 



A. Rotation frequency versus hole radius Q over 

Eliminating Rq from the il (Ro) relation for free spirals, 
Eq. ( 49 I by using Eq. ( 54 1 , we obtain 



R+n = R+n 1 - 




41^2 



(i?+0- 1)2 



(83) 



Because i?+ is the given hole radius, (83) already 
represents the desired result. Note that in contrast to 
free spiral waves, for a pinned spiral wave the rotation 
frequency does not depend on the parameter B that 
characterizes the strength of the interaction between 
wave front and back and measures the excitability of 
the medium. Within our approach, any dependence of 
the rotation frequency on B or other kinetic parameters 
enters through the dispersion relation for periodic pulse 
trains. 

As before, for small e an explicit expression can be de- 
rived perturbatively 



oj (r+) = 



r+ 



2ci/4 

7/4 ' 



3/4 



+ O . (84) 
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Figure 6: Rotation frequency SI versus hole radius R+ for 
spiral waves pinned to a Neumann hole. The analytical result 
(black solid line, Eq. ( |83[ )) is compared to Keener's result 
(green dot-dashed line, the result by Hakim and Karma 
(blue dotted line, |14|'). and numerical solution of the linear 
eikonal equation (red dashed line). 
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Eq. (84 1 can be compared with Keener's result ^ 
c ('4cr+ + e - ye(8c7>Te)) 

= , 

which after expanding for small e gives 

cl/2 



WK (r+) = 



.1/2 



/2r: 



3/2 



0(e) 



In rescaled form, Keener's result reads 

1 + 4i?+ - yr+ 



{R+) = 



8B., 



4i?+ (1 + 



(85) 



(86) 



(87) 



Because in general the difference between the core ra- 
dius tq and the radius of the corresponding Neumann 
hole r_|_ is small, Keener's result can be applied to free spi- 
ral waves. This was done successfully by Winfree, com- 
pare The result obtained by Hakim and Karma 
can be modified for pinned spirals according to 



^HK {r+) = h 



5/3 



.2/3, 



where here ai = —1.01879 denotes the global maximum 
of the Airy function Ai(a;). The rescaled form of this 
expression 



^HK(i?.+ ) = ^+2V3_^, 



R 



R 



(89) 



together with our result for Vl{R+) and Keener's result 
f^K {R+); is compared in Fig. plwith numerical simula- 
tions of the rescaled eikonal equations Eqs. (31 1, (32 1. 



Figure 7: Shape of spiral waves pinned to a Neumann hole. 
Analytical result (black solid and dashed lines) and numerical 
results (red dotted and dot-dashed lines) are plotted for dif- 
ferent hole radii and given excitability: a) i?+ = 2, B = 0.4, 
b) i?+ = 4, B = 0.3, c) i?+ = 3, B = 0.7. 



All three analytical results agree well with the numer- 
ically obtained curve in the limit of large hole radius 



i?4 



cxD while for small core radii i?, i 



partially 



marked deviations appear. Merely our analytical approx- 
imation produces a finite and different from zero rotation 
frequency for vanishing hole radius. Note the different ex- 
ponents in leading order of e in the expansions Eq. ( 84 1 , 
Eq. ([86|, and Eq. (ISSl. 



B. Approximations for the wave shape Bans {R) 

From the Neumann boundary condition at i? = i?,+ , we 
have R- = i?+ and it follows 



*^„s [R) - *ans {R) ■ 

For the shape of the front interface we obtain 



(90) 



4Q2 

0+ . {R) = Qt-n,lR^ + R\[ 1 



+ 



Rj^Q. arccos \l 



1 - 



(fl,+ 0-l)2 



1 



(91) 



where denotes an arbitrary initial phase. The relation 
between hole radius and rotation frequency is given 
by Eq. (83 I. Furthermore, from asymptotic solutions of 



the rescaled eikonal equations Eqs. (31 1, (32 1 at i? = 



11 



2n 

we find 0+ (i?+) — {R+) = -^r, from which follows 

B 



©Ins (R) = 0a"„s (R) + 



2n 

l3' 



(92) 



We compare the analytical predictions according to Eq. 



(91 1 with numerical solutions of the kinematic equations 
in Fig. [Tj Close to the hole, the agreement between the- 
ory and numerics is good. With larger R values, initially 
minor deviations arise which grow according to 



e± (i?)-e±,(i?) = 7i? + o(i?), R 



oo, 



(93) 



where 7 is some nonzero constant. The reason for this 
discrepancy is that, first, the analytical relation 51 (-R+) is 
only approximately valid, and, second, the correct lead- 
ing order asymptotics for the shape functions is given by 
e± (i?) -VtR + O (log (i?)) , i? ^ 00. 



VI. CONCLUSIONS 



Based on a new non-perturbative ansatz Eq. ( 40 1 , we 



have presented analytical approximations for spiral wave 
solutions to the linear eikonal equation. The approx- 
imate analytical solution (r) for a rigidly rotating 
spiral wave displays the correct leading order asymptotic 
expansion of the unknown exact solution 6^ (r) close to 
the core (r — > tq), at the radius of the effective Neumann 
hole (r — >■ r+, 0^' (r)) and far from the core (r — >■ 00). In 
addition, approximate and exact solution have the same 
asymptotic expansion in the inflection point located on 
the wave back (r — > ri, k (ri) = 0). The asymptotically 
correct treatment of these four essential regions of a spi- 
ral wave is crucial for our ansatz. 

For the front interface our ansatz works quite good. In 
particular, the derived dependence between the rotation 
frequency Vt and the core radius Rq for rigidly rotating 
spiral waves, Eq. (49 1, agrees well with numerical so- 
lutions of the linear eikonal equation. In fact, for large 
and intermediate core radii the agreement between the- 
ory and numerics is very close, moreover, even for small 
core radius our ansatz produces a more than acceptable 
match with the numerical results. An equally good ana- 
lytical approximation for f2 (i?o)i which is globally valid 



for all spiral core radii i?o, to our knowledge, does not 
exist. 

Using singular perturbation theory it is impossible to de- 
rive a globally valid approximation for lo (ro) starting 
from Eq. ( [48| . The reason is that the small parame- 
ter e by introducing rescaled quantities drops out of Eq. 
). This scaling gives a dominant balance. There is no 
ing giving a dominant balance which leads from Eq. 
(|48| to an equation which is simpler than Eq. (49 1. In 



(49 



other words, sooner or later we will inevitably be faced 



with Eq. ( 48 1 in order to obtain a globally valid solution 
for u) (ro). 

We believe the same to be true for the full free boundary 
approach based on the linear eikonal equations Eqs. (21 1, 



( 22 1 because they can be transformed to the rescaled 



eikonal equations Eqs. (31 1, (32 1 that do not contain 



the small parameter e any more. 
For the shape of the back interface and the dependence of 
the rotation frequency Q. (B) on the excitability parame- 
ter B, the results produced by the ansatz are much less 
satisfactory. The reasons for the discrepancy between an- 
alytical prediction and numerical results are discussed in 
sections HVB] and HVCl 

The relatively simple i?-dependence for the wave shape 



Eq. ( 69 ) could be the starting point for stability analysis 



and further analytical investigations. 
In principle, one can improve the ansatz Eq. (41 1 by 
including terms involving additional constants. Higher 
order asymptotics at tq, r^, ri and for r — >■ (X) can be 
taken into account to determine these constants. How- 
ever, it is difficult to find terms which show the correct 
asymptotic behavior at one point without simultaneously 
destroying the correct asymptotics at other points. 
Rotation frequency and core radius are accessible in ex- 
periments with the Belousov-Zhabotinsky reaction, for 
example. The relation for Vl{Rq), Eq. (49 1, supple- 



mented with the dispersion relation for one-dimensional 
pulse trains, can be checked in experiments. 
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